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ABSTRACT 

Gamma-ray emission from pulsars is thought to arise from accelerating regions in 
pulsar's outer magnetosphere. The shape of the light curves is thus sensitive to the 
details of the magnetic geometry of the magnetosphere. In this work, we show the first 
calculations of light curves from the more realistic force-free field under the framework 
of conventional emission models. We compare the properties of gamma-ray emission 
between the commonly used vacuum dipole magnetic field and the new force-free field. 
We discuss the role of the polar cap shape and aberration effect on the appearance of 
the light curves as well as the formation of caustics on the sky map. With the force- free 
field, the double-peak pulse profile is best reproduced if the emission zone lies in a 
thin layer just outside the current sheet, and the peaks are mainly contributed from 
regions near the light cylinder. The conventional two-pole caustic model can produce 
up to four peaks in general, while the conventional outer-gap model can normally 
produce only one peak. These results will be useful for interpreting Fermi telescope 
observations. 
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1 INTRODUCTION 

Pulsars are rotating neutron stars (NS) with strong mag- 
netic field. Seven of them were detected as gamma-ray pul- 
sars with EGRET (Thompson 2004). The light curves of 
these gamma-ray pulsars are typically double-peaked, and 
have substantial off-peak emission. Theoretical models have 
been developed in order to explain the nature of pulsar 
gamma-ray emission, namely, the polar-cap (PC) model 
(Daugherty & Harding 1982; Daugherty & Harding 1996), 
the slot-gap (SG, or two-pole caustic, TPC for short) model 
(?; Arons 1983; Dyks & Rudak 2003; Dyks et al. 2004), and 
the outer-gap (OG) model (Cheng et al. 1986; Romani & 
Yadigaroglu 1995; Cheng et al. 2000). In these models, par- 
ticles are accelerated to ultra-relativistic energies in "gap" 
regions where strong electric fields are developed due to the 
deficit of charge. The gamma ray emission comes from curva- 
ture and inverse Compton (IC) radiation of these energetic 
particles. Based on vacuum magnetic field geometries, SG 
and OG models are more favored since they can reasonably 
well reproduce the double-peak light curves. 

However, pulsar magnetosphere is filled with plasma 
(Goldreich & Julian 1969). The plasma is essentially force- 
free (FF), with pEi + jxB/c = 0, where p and j are charge 
and current density. In the presence of plasma, the magneto- 
sphere consists of open and closed field line regions. Poloidal 
current flows out along open field lines and induces toroidal 
magnetic field, which dominates beyond the light cylinder 



(LC for short, R LC = c/fi). Recently, Spitkovsky (2006) ob- 
tained the full 3D magnetospheric structure using FF sim- 
ulations. The FF field geometry differs substantially from 
the vacuum field geometry (i.e., the retarded dipole field), 
which is commonly used in light curve calculations, espe- 
cially around the LC. The light curve is very sensitive to the 
geometry of the emission zones, thus to the field geometry 
itself. Therefore, it is important to revisit pulsar high-energy 
emission models using the more realistic FF field geometry. 

This paper is organized as follows. In section 2 we show 
the light curves from the vacuum field geometry, pointing 
out that previously calculated vacuum light curves should 
be modified when aberration effect is correctly treated. In 
section 3 we present the light curves from the FF field using 
the two-pole caustic model, and mention the applications to 
the outer-gap model. 



2 RELATIVISTIC EFFECTS AND VACUUM 
FIELD RESULT CORRECTIONS 

In order to calculate the light curve, we need to find the 
emitting region in the magnetosphere. In the TPC model, 
the emission zone is assumed to be along the last open field 
lines (LOFLs) starting from the polar cap. In the OG model, 
the radiation region is assumed to be beyond the null-charge 
surface and along open field lines. Technically, given a mag- 
netic field geometry (either vacuum or FF field), field lines 
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are considered open if they can cross the LC in the lab (ob- 
server's) frame (LF). We trace magnetic field lines in the 
LF and find the LOFLs. The polar cap is the region on the 
NS surface where open field lines originate. We calculate the 
magnetic colatitude of polar cap rim 0™ at fixed magnetic 
azimuth m (the subscript "m" means w.r.t. magnetic axis 
rather than rotation axis), and define the open volume coor- 
dinate as r ov = 6> m /6>™. In the TPC model (Dyks & Rudak 
2003; Dyks et al. 2004) emission comes from a thin layer 
centered at LOFL (r ov = 1) extended from the NS surface 
up to a certain radius (R ~ 0.75i?LC,f ~ I.CLRlc where r 
is distance to the NS center and R is cylindrical radius). In 
the OG model (Cheng et al. 2000) the emission comes from 
a thicker layer centered at r ov ~ 0.90 extended from the null 
charge surface (where B z = 0) to the LC. 

For curvature/IC radiation, the radiation direction 
should be along the direction of particle motion. An ultra- 
relativistic particle moves at the speed of light, and its veloc- 
ity can be decomposed into corotation velocity (i.e., drift ve- 
locity) plus a component along the direction of the magnetic 
field. Thus, the radiation direction ff in the LF is determined 
by 

if=fB+0, (1) 

where /3 — Q x r/c is the normalized corotation velocity, 
B is the magnetic field in the LF, and / is a coefficient 
determined by the requirement that \ff\ = 1. 

This approach is different from earlier works on OG and 
TPC models with vacuum field geometries [e.g. Romani & 
Yadigaroglu (1995), Cheng et al. (2000), Dyks et al.(2004)]. 
In previous works, it was implicitly assumed that the vac- 
uum field, which is the solution of vacuum Maxwell equa- 
tions in the LF, is valid in the instantaneously corotating 
frame (ICF). In other words, in eqn.(l), B was improperly 
replaced by B', where B' is obtained by the Lorentz transfor- 
mation of B from ICF to LF [see Bai k Spitkovsky (2009a) 
for a more detailed discussion]. Therefore, results from these 
models require revision. 

Photons emitted from the emission zones are collected 
in the sky map (0,£ o bs), where <j> is the phase of rotation 
with corrections for photon travel time (i.e., time delay ef- 
fect), and £obs is observer's viewing angle. We assume con- 
stant emissivity along field line. The light curve seen by the 
observer is then obtained by cutting the sky map at a spec- 
ified viewing angle £ bs- 

In Fig.l we plot the sky map and light curve from 
TPC model for an oblique rotator with the inclination an- 
gle a — 60° using vacuum field geometry. The upper plot is 
a reproduction of Dyks et al. (2004) 's result, where the re- 
tarded vacuum field was implicitly assumed to be in the ICF. 
There are two strong caustics associated with the two poles. 
They are formed at modest distance from the NS (roughly 
0.3 — 0.6i?Lc)- In the light curve, two sharp peaks are thus 
produced. In the bottom plot we correct the aberration with 
eqn.(l). We find that the caustics in the upper plot are now 
much widened and weakened. There are still two peaks, but 
much wider, and a substantial contribution to the peaks 
comes from the overlap of emission from both poles, rather 
than from the caustics. This result also holds for other in- 
clination angles. We conclude that the TPC model has diffi- 
culties in producing sharp peaks with corrected aberration. 
We also note that the correction of abberation formula also 
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Figure 1. Sky map (left) and light curve (right) from two-pole 
caustic model using vacuum field. Inclination angle a = 60°, 
viewing angle £ b s — 80°. Up: with inappropriate treatment of 
aberration; down: with the correct aberration [eqn.(l)]. formula. 



affect the conventional outer-gap models using the vacuum 
field (Bai & Spitkovsky 2009a). 



3 RESULTS FROM FORCE-FREE FIELD 
GEOMETRY 

The advent of FF field (Spitkovsky 2006) makes it possi- 
ble to test theoretical models with more realistic magneto- 
spheric geometries. The FF magnet osphere has E • B = 
everywhere and thus has no intrinsic particle acceleration. 
As long as pulsar's emission power is much smaller than the 
total spin down energy loss rate (which is almost always the 
case), the FF field provides a reliable field structure in the 
magnet osphere. The pulsar's radiation comes from a small 
region of magnetosphere (i.e., the emission zone) where non- 
ideal deviations from the FF condition makes particle accel- 
eration possible. The current sheet is a key feature of FF 
field. Beyond the LC, it separates field lines from the two 
poles (Contopoulos et al. 1999; Bai & Spitkovsky 2009b), 
and it is also a likely place for dissipation and resistive heat- 
ing (Gruzinov 2007). Consequently, the pulsar's gamma-ray 
emission may be connected to the current sheet. 

The sky map is determined by the field geometry as well 
as the shape of the polar cap. Although the FF field appears 
similar to the vacuum field near the NS surface, the shape 
of the polar cap is different, because the polar cap shape is 
sensitive to the field structure near the LC. The polar cap 
shape in the FF field is more circular and larger than in the 
vacuum field (Bai & Spitkovsky 2009b), and it will lead to 
fundamental differences in the sky maps, as we will illustrate 
here. 

In Fig. 2 we show the sky map of the conventional TPC 
model using the FF field which is calculated with r ov = 1.0. 
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Figure 2. Sky map from conventional two-pole caustic model 
r ov = 1.0 using force- free field. Inclination angle a = 60°, viewing 
angle Cobs = 80°. 




Figure 3. Sky map (left) and light curve (right) from the annular 
gap model (a modified two-pole caustic model) with r ov =0.90 
(up) and the conventional outer-gap model with r ov =0.90 (bot- 
tom) using force-free field. Inclination angle a = 60°, viewing 
angle £ obs = 80°. 



The emission zones are assumed to be extended from the 
polar cap to the LC along current sheet. Depending on ob- 
server's viewing angle, the sky map can produce light curves 
with up to four widely separated peaks, which are inconsis- 
tent with the current observations. We have experimented 
with a number of other possible solutions, and find that the 
double-peaked pulse profile can be reproduced if we allow 
the emission zones to be located in a layer just outside the 
current sheet. In Fig. 3 (upper plot) we show the sky map 
and light curve from the emission zone Tov = 0.90 under 
this "annular gap" model (or modified TPC model). There 
are now two strong caustics associated with both poles. The 
caustics here are formed around the LC. The double-peak 
pattern is very robust as we vary the inclination angle a. 

The choice of the open volume coordinate r ov = 0.90 is 
similar to that used in the OG models. However, we empha- 
size here that this is not an OG model. For an OG model, 
only field lines that cross the null- charge surface can con- 



tribute (Takata et al. 2008). Here we have considered emis- 
sions from all field lines with a certain r ov - It turns out that 
at ^bs = 60°, only a small fraction of the field lines can 
cross the null-charge surface, and only one peak survives 
(see bottom plot of Fig. 3). 



4 CONCLUSION 

Existing models for pulsar's gamma-ray emission (i.e., slot- 
gap and outer-gap models) are all based on the vacuum field 
geometries, which carry large uncertainties. Hence, the light 
curves predicted from these models are questionable. In this 
paper we revisit pulsar's gamma-ray emission using the more 
realistic force- free field configurations. 

We first point out the inconsistent treatment of the 
aberration effect in previous works and compare them with 
aberration corrected results. The new results indicate that 
sharp peaks are difficult to be produced using vacuum field 
geometry in the conventional two-pole caustic model. 

Using the force- free field, we find that the conventional 
slot-gap (two-pole caustic) model tends to produce four 
peaks, while the conventional outer-gap model is able to 
provide just one peak, inconsistent with current observa- 
tions. An "annular gap" model, or modified two-pole caustic 
model, where emissions are assumed to originate from a thin 
layer just outside the entire current sheet, is able to repro- 
duce the basic observed features of gamma-ray pulsars. The 
physics behind remains to be understood. 

The launch of the Fermi Gamma-ray Space Telescope 
will significantly expand the sample of gamma-ray pulsars 
and will provide very accurate pulse profiles. Our results will 
be useful for the interpretation of Fermi telescope observa- 
tions and for probing the physics of pulsar magnet ospheres. 
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